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Abstract 

We analyze the connections between the mathematical theory of 
knots and quantum physics by addressing a number of algorithmic 
questions related to both knots and braid groups. 

Knots can be distinguished by means of 'knot invariants', among 
which the Jones polynomial plays a prominent role, since it can be 
associated with observables in topological quantum field theory. 

Although the problem of computing the Jones polynomial is in- 
tractable in the framework of classical complexity theory, it has been 
recently recognized that a quantum computer is capable of approxi- 
mating it in an efficient way. The quantum algorithms discussed here 
represent a breakthrough for quantum computation, since approximat- 
ing the Jones polynomial is actually a 'universal problem', namely the 
hardest problem that a quantum computer can efficiently handle. 

Key words: knot theory; Jones polynomial; braid group representations; spin- 
network simulator; quantum computation. 
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1 Introduction 



Knots and braids, beside being fascinating mathematical objects, are en- 
coded in the foundations of a number of physical theories, either as concrete 
realizations of natural systems or as conceptual tools. The atomic model 
based on knot theory, proposed in the nineteenth century, is a well known 
(although wrong) example of interaction between the experimental analysis 
of reality, on the one hand, and the need for classifying the observed struc- 
tures in well defined mathematical categories, on the other. A much more 
recent example is provided by 'knotted' configurations of DNA strands. 
Since knots are collections of 'knotted' circles and braids are 'weaving' pat- 
terns, the term 'knitting' in the title of this paper seems the more appro- 
priate to describe collectively these structures. Moreover, knots and braids 
are indeed closely interconnected since from a braid we can get a multi- 
component knot (or a 'link') by simply tying up its free endpoints, while 
from a knot drawn in a plane we can select portions which look like 'over' 
and 'under '-crossing strings (c/r. Fig. ^and Fig. |1] below). 

On the conceptual side, it was in the late 1980 that knot theory was rec- 
ognized to have a deep, unexpected interaction with quantum field theory 
pp. In earlier periods of the history of science, geometry and physics inter- 
acted very strongly at the 'classical' level (as in Einstein's General Relativity 
theory), but the main feature of this new, 'quantum' connection is the fact 
that geometry is involved in a global and not purely local way, i.e. 'topolog- 
ical' features do matter. To illustrate this point, consider a knot embedded 
in the three-dimensional space. What really matters is the way in which it 
is 'knotted', while the length of the string and the bending of the various 
portions of the string itself can be changed at will (without cutting and 
gluing back the endpoints). Over the years mathematicians have proposed 
a number of 'knot invariants' aimed to classify systematically all possible 
knots. As will be illustrated in sections 2 and 3, most of these 'invariants' 
(depending only on the topological features of the knot) are polynomial ex- 
pressions (in one or two variables) with coefficients in the relative integers. 
It was Vaughan Jones in |21 who discovered the most famous polynomial 
invariant, the Jones polynomial, which connects knots with quantum field 
theory and plays a prominent role in the present paper. In the seminal work 
by Edward Witten [1^, the Jones polynomial was actually recognized to be 
associated with the vacuum expectation value of a 'Wilson loop operator' 
in a particular type of three-dimensional quantum field theory (the non- 
abelian Chern-Simons theory with gauge group SU{2)). 
Braids appear naturally in this context too, since we can always 'present' 
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a knot as the closure of a braid. Moreover, braids enrich the purely topo- 
logical nature of the theory since the set of crossings of any braid can be 
endowed with a group structure. The Artin braid group on n strands (to 
be defined in section 2) encodes all topological information about 'over' and 
'under' crossings into an algebraic setting, opening the possibility of describ- 
ing polynomial invariants of knots in terms of 'representations' of this group. 
Thus, the Jones polynomial can be interpreted not only as the 'trace' of a 
(suitably chosen) matrix representation of the braid group [2] , but also as the 
operatorial trace of an observable in a unitary quantum field theory In 
section 3 below we shall provide the (minimum) mathematical background 
needed for understanding this crucial issue, and realizing 'quantum knitting'. 

The search for new algorithmic problems and techniques which should 
improve 'quantum' with respect to classical computation is getting more and 
more challenging. Most quantum algorithms run on the standard quantum 
circuit model |3j, and are designed to solve problems which are essentially 
number theoretic (such as the Shor algorithm). However, other types of 
problems, typically classified in the field of enumerative combinatorics, and 
ubiquitous in many areas of mathematics and physics, share the feature to be 
'intractable' in the framework of classical information theory. In this latter 
perspective, the first section of the present paper deals with a list of algorith- 
mic questions in knot theory and in the theory of finitely presented groups, 
focusing in particular on the braid group and its algebraic structure encoded 
into the Yang-Baxter relation. These examples are interesting not only for 
mathematicians and computer scientists, but are addressed by physicists in 
the study of both exact solvable models in classical statistical mechanics 
[1] (where Yang-Baxter equation reflects the integrability properies of such 
models) and quantum field theories (such as topological field theories men- 
tioned above and conformal field theories [Hj). 

In section 3 we give the essential mathematical apparatus underlying the 
construction of link polynomials. This will set the stage for the last section, 
where algorithmic questions concerning the calculation of such topological 
invariants are raised and discussed. Indeed, in the late 2005 Dorit Aharonov, 
Vaughan Jones and Zeph Landau 6^ found an efficient quantum algorithm 
for approximating the Jones polynomial (a difficult problem to be addressed 
in the classical context), thus giving a fundamental improvement of pre- 
vious claims on the computational complexity of such problem (Jj. Two 
papers of ours [HI IH] provided essentially similar results for a wider class 
of link invariants, and the paper by Wocjan and Yard ^Uj addressed and 
extended the algebraic and information-theoretical background used in |^. 
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The quantum algorithm in question (reviewed and discussed in section 3) 
represents a breakthrough for quantum computation because the problem 
of approximating the Jones polynomial has been recognized to be a 'uni- 
versal problem', namely the hardest problem that a quantum computer can 
efficiently handle. The latter remark sounds reminiscent of 'quantum knit- 
ting': we can really knit knots and braids by means of 'quantum' computing 
machines. 

2 Knots, braids and related classical 
algorithmic problems 

A knot K is defined as a continuous embedding of the circle (the 1- 
dimensional sphere) into the euclidean 3-space or, equivalently, into the 
3-sphere = U {oo}. A link L is the embedding of the disjoint union 
of M circles, ^^=1 {S^)m into or 5*^, namely a finite collection of knots 
referred to as the components of L and denoted by {Lm}m=i,2,...,M- Since 
each circle can be naturally endowed with an orientation, we can introduce 
naturally oriented knots (links). 

Referring for simplicity to the unoriented case, two knots Ki and K2 are 
said to be equivalent, Ki ~ if and only if they are (ambient) isotopic. 
An isotopy can be thought of as a continuous deformation of the shape of, 
say, K2 C which makes K2 identical to Ki without cutting and gluing 
back the 'closed string' K2. 

The planar diagram, or simply the diagram, of a knot K is the projection 
of -ftT on a plane C M'^, in such a way that no point belongs to the 
projection of three segments, namely the singular points in the diagram are 
only transverse double points. Such a projection, together with 'over' and 
'under' information at the crossing points -depicted in figures by breaks in 
the under-passing segments- is denoted by D{K); a link diagram D{L) is 
defined similarly. In what follows we shall sometimes identify the symbols K 
[L] with D{K) [D{L)], although we can obviously associate with a same knot 
(link) an infinity of planar diagrams. Examples of diagrams are depicted in 
Fig. □ 

After these preliminary definitions, let us state the first of our question 
on algorithmic complexity, actually the fundamental, still unsolved problem 
in knot theory. 

Problem . Give an effective algorithm for establishing when two 
knots or links are equivalent. 
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Figure 1: Planar diagrams of the trefoil knot (top) and Borromean link 
(bottom). 

The number of crossings of a knot (diagram) is clearly a good indicator of the 
'complexity' of the knot. Indeed, Tait in late 1800 initiated a program aimed 
to classifying sistematically knots in terms of the number of crossings (see 
jllj , |12[ , IJl] for exhaustive accounts on knot theory and for references 
to both older papers and knots tables). 

Since a knot with crossing number n = c{K) can be represented by planar 
diagrams with crossing numbers c{D{K)) for each c{D{K)) > k, the first 
issue to be addressed is the search for procedures aimed to simplify as much 
as possible the diagrams of a knot K to get a D'{K) with c{D'{K)) = k, 
the 'minimum' crossing number. Reidemeister's theorem (see |12) ) gives the 
answer to this basic question. 

Equivalence of knots (Reidemeister moves). Given any 
pair of planar diagrams D,D' of the same knot (or link), there 
exists a finite sequence of diagrams 

D = Di ^ D2 ^ > Dk = D' (1) 

such that any -Dj+i in the sequence is obtained from Di by means 
of the three Reidemeister moves (I, II, III) depicted in Fig. [2 

The procedure determined by Reidemeister's theorem applies to subsets 
of link diagrams localized inside disks belonging to the plane where the 
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Figure 2: The three Reidemeister moves acting on local configurations in 
link diagrams. 

diagram lives, and can be suitably generalized to handle diagrams of oriented 
links. However, notwithstanding the recursively numerable character of the 
implementation of the Reidmeister moves with respect to the intractability of 
the notion of ambient isotopy, such moves cannot be formalized into effective 
algorithms, basically because the above definition is purely topological. As 
we shall see, transformations on link diagrams can be consistently translated 
into an algebraic setting by exploiting their deep connection with braid 
groups. In the new setting the 'moves' on link diagrams will be reformulated 
in terms of algebraic operations, see Markov's theorem below. 

Before addressing this issue, which will generate several interesting algo- 
rithmic problems, let us point out that in geometry (as well as in physical 
field theories whose dynamical variables have a geometric nature) an impor- 
tant role is played by 'topological invariants'. To illustrate this point and 
clarify its connection with the systematic classification of geometric struc- 
tures, consider the case of surfaces. In particular, any smooth, closed and 
oriented surface S has a well defined topological structure, completely de- 
termined by its Euler number x{<S) = 2 — 2g, where g is the number of 
'handles' of S. Thus we have the 2-sphere S'^ for g = (no handles), the 
2-torus (the surface of a doughnut) for (7 = 1, and more complicated, many- 
handles surfaces as g grows. The existence of such a 'complete' invariant 
in 2-dimensional geometry is quite an exception: higher dimensional geo- 
metrical objects (smooth D-manifolds) do not share this feature, and the 
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associated classification problems are still open (except for some restricted 
classes). 

Coming back to knot theory, we define a link invariant through a map 



L 



f{L) 



(2) 



where the quantity f{L) depends only on the type of the link, namely takes 
different values on inequivalent links. Switching to link diagrams, we keep 
on using the same notation as in 0, but now it is sufficient to verify that 
f{L) (= f{D{L))) does not change under applications of the Reidemeister 
moves I, II, III (Fig. 

We have already met a numerical invariant, namely the (minimum) cross- 
ing number k. It is a natural number which takes the value for the trivial 
knot represented as an unknotted circle. Other invariants taking values in 
Z (relative intergers) can be defined for oriented link diagrams, where each 
crossing is marked by ±1 according to some fixed convention. For instance, 
the writhe w{D[L)) of a diagram D of an oriented link L is the summation 
of the signs of the crossings of namely 



where the sum runs over the crossing points {p} and ep = +1 if the (di- 
rected) knot path shows an overpass at the crossing point p, tp = —1 for an 
underpass. Note however that both the crossing number and the writhe do 
change under Reidmeister move of type I, but are invariant under the moves 
II and III: this property defines a restricted kind of isotopy, commonly re- 
ferred to as regular isotopy. The concept of regular isotopy is very useful 
because, by eliminating the move I, we do not really lose any information 
about the topology of the link. Moreover, the evaluation of crossing num- 
bers and writhes can be carried out efficiently by a simple inspection of the 
diagrams. 

Over the years, mathematicians have provided a number of knot in- 
variants, by resorting to topological, combinatorial and algebraic methods. 
Nevertheless, we do not have yet a complete invariant (nor a complete set 
of invariants) able to characterize the topological type of each knot and to 
distinguish among all possible inequivalent knots, c/r. Problem 0. 
As a matter of fact, the most effective invariants have an algebraic origin, 
being closely related to the braid group and its representation theory. As we 
shall see below, it is straightforward to obtain a knot (link) out of a 'braid'. 
The inverse process is governed by Alexander's theorem (see jTH], section 2). 




(3) 



p 



7 



Braids from links. Every knot or link in 5''^ = R'^ U {00} can 
be represented as a closed braid, although not in a unique way. 

The Artin braid group B„, whose elements are (open) braids /5, is a finitely 
presented group on n 'standard' generators {ai, cj2, . . . , ct„_i} plus the iden- 
tity element e, which satisfy the relations 

o-jcjj =ajai (i, j = 1,2, . . . ,n - 1) if \i - j\ > 1 
(Tj o-j+i cJi = (Tj+i CTj cJi+i (z = 1,2, . . . ,n - 2). (4) 

This group acts naturally on topological sets of n disjoint strands with fixed 
endpoints - running downward and labeled from left to right - in the sense 
that each generator fjj corresponds to a crossing of two contiguous strands 
labeled by i and (i + 1), respectively: if ai stands for the crossing of the i-th 
strand over the (i + l)-th one, then a~ represents the inverse operation 
with (Tj a^^ = (T~^(Ti = e, see Fig. |21 An element of the braid group can be 
thought of as a 'word', such as for instance /? = ^^^^2 a^^a2(y\ ct^"^(Ti(T^^ 
G B4; the length \(3\ of the word (3 is the number of its letters, where by a 
'letter' we mean one of the generators or its inverse element. 
By a slight change of notation, denote by the over-crossing operation 
acting on two strands the endpoints of which are labelled by i and j. Then 
the second relation in (jlj can be recasted into the form 

R12 Ri3 R23 = R23 Ri3 R12 (5) 

and represented pictorially as in Fig. |31 where operations are ordered down- 
ward. Note that this picture can be viewed as a portion of an n-strands 
configuration (and thus {1,2,3} may actually represent labels attached to 
any triad of contiguous strands) since the first relation in (0]) ensures that 
other types of crossings cannot happen at all. The relation (0) is referred 
to as the algebraic Yang-Baxter relation and characterizes the structure of 
solvable models in statistical mechanics jU |3] . 

As mentioned before, it is straightforward to get a link out of a braid: 
we have simply to 'close up' the ends of an open braid (3 to get a closed braid 
(3 that reproduces the diagram of some link L. Formally 

/3 /3 ^ L. (6) 

Notice however that this operation can be performed in two ways, denoted 
by (the standard closure) and /3p' (the plat closure), respectively. In 
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Figure 3: The generator cjj and its inverse (top). The algebraic Yang- 
Baxter equation (bottom). 

Fig. the two admissible closures of a same open braid are shown, and 
moreover the content of Alexander's theorem is made manifest, since both 
these closed braids can be seen as deformations (isotopy transformations) of 
the planar diagram of the trefoil knot depicted in Fig. ^ 

As already pointed out, Alexander's theorem does not establish a one- 
to-one correspondence between links and braids. For instance, given a closed 
braid (3 = L with (3 € B„, any other braid obtained from (3 by conjugation, 
namely /?' = a/3a~^ (for some a G B„) has a closure (3' which reproduces 
the same link L. Thus the following question can be naturally rised. 

Problem lA . Is it always possible to transform efficiently a given 
knot or link into a closed braid? 

The answer is affirmative, since there exists a classical algorithm which 
performs the reduction in a number of steps which is bounded from above 
by a polynomial function of the braid index (|13|. section 2 and original 
references therein) , where the braid index of a braid or closed braid is simply 
the number of its strands. 

Taken for granted the above result about the efficiency of the reduction 
of any link diagram to a closed braid, we can now exploit the algebraic 
properties of braid groups. For what concerns in particular the issue of 
equivalence, Reidemeister's theorem can be recasted into Markov's theorem. 
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Figure 4: The two types of closures of a braid, namely the plat closure (left) 
and the standard closure (right), both representing the same trefoil knot. 

The following statement of this theorem refers to the case of open braids, 
which captures the crucial features of the construction, while the version 
involving closed braids can be found in |13j . section 2. 

Equivalence of braids (Markov moves). Two braids are 
equivalent if they differ by a finite sequence of Markov 'moves' 
of the following two types, together with their inverse moves: 

i) change a braid (3 G B„ to a conjugate element in the same 
group, (3 — > af3a~^, with a G B„; 

ii) change /? G B„ to i„(/?)(T^^, where i„ : B„ B„_|_i is the 
natural inclusion obtained by disregarding the (n + l)th strand 
and cr„,cr^^ G B„. 

The next question arises in connection with the search for the most 'eco- 
nomical' representation of a knot diagram as a closed braid. The minimum 
braid index of a link L is the minimum number n for which there exists a 
braid (3 G B„ whose closure /3 represents L. 

Problem IB . Does there exist an (efficient) algorithm to select, 
among all the diagrams of a given link L, the diagram with the 
minimum braid index? 

At present no explicit algorithm for addressing this problem is known, so 
that its computational complexity class cannot be even evaluated (see 
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section 2 and 4 for more details). 

Coming to algorithmic problems characterized in purely algebraic terms, 
let us recall that braid groups belong to the class of finitely presented groups. 
Such groups are defined by means of a finite sets of generators together 
with relations among the generators and can have a finite order -as for 
the classical point groups of cristallography- or not -as happens for the 
braid group- (see and older references therein). It was Max Dehn who 
stated the three 'fundamental problems' (listed below) concerning a group 
G presented in terms of generators, denoted by a,b,c,..., and relations 
P,Q, R, . . ., namely 

G = {a,b,c,... ■,P,Q,R,...) . (7) 

Any element of G can be written (in multiplicative notation) as a word W in 
the alphabeth given by the generators and their inverse elements, as already 
done for the braid group. Note that the relations P, Q,R, ... in ^ represent 
the minimal set of words in the generators, equivalent to the identity element 
e ('minimal' meaning that any other word equivalent to the identity can be 
reduced -by the use of the relations in the set- to the union of words in the 
same set). 

A simple example of relation for a group G on two generators, a, b, is ab = ba 
or, equivalently, P = ab{ab)~^ = e. Of course this means that the group is 
commutative, and looking back at the presentation of the braid group given 
in @ we see that the generators fjj and Uj do commute when |i — i| > 1, 
but B„ is not a commutative group since (TjiTi+i ^ crj+ifTj. 
Given the presentation (O of the group G, Dehn's problems are formulated 
as follows. 

2A. The word problem. For an arbitrary word W in the genera- 
tors, decide whether or not W defines the identity element in G. 
Equivalently: given two words W , W , decide whether W = W . 

2B. The conjugacy problem. For two arbitrary words Wi and W2 
in the generators, decide whether or not they are conjugate to 
each other. In a sharper form: find explicitly an element W for 
which W2 = W'WiiW')-^. 

20. The isomorphism problem. For an arbitrary group G' de- 
fined by another presentation G' = {a', b' , c', . . . ; P', Q', R', . . .), 
decide whether or not G and G' are isomorphic. 
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Except for the second issue in Problem 2B, we are in the presence of deci- 
sion problems, namely problems that can be addressed by means of classical, 
deterministic algorithms (running on a Turing machine) designed to answer 
'yes' or 'no' to each of the above questions. Recall also that the time com- 
plexity function f_4 of an algorithm A is defined in terms of the size (length) s 
of the input. In particular, an algorithm associated with a decision problem 
belongs to the complexity class P if f_4 (s) is a polynomial function of s and 
to the class NP is any guess on the answer can be checked in polynomial 
time. Algorithmic problems endowed with complexity functions of exponen- 
tial type are to be considered as intractable in the framework of classical 
information theory. 

For what concerns the list above in the case of a generic group G, note 
first of all that a solution of the complete conjugacy problem 2B (belonging 
to a certain complexity class) would imply a solution to Problem 2A (in the 
same class) since it would be sufficient to set W' = e in the expression W2 = 
W'Wi{W')~^ . It also clear that the most difficult problem is the last one, 
which requires a 'global' inspection of the algebraic structures (generators + 
relations) of the groups under examination. As for the braid group, we leave 
aside this problem and refer the reader to (section 1) for the definition 
of a presentation in terms of generators and relations alternative to the 
standard ones collected in Q. 

The known results about the word and the conjugacy problems for the 
braid group B„ are briefly summarized below (see section 5 and original 
references therein). 

2A. The solution to the word problem is polynomial, with a 
complexity function of the order 



where |/3| is the length of the initial representative of the braid 
P and n is the braid index. 

Surprisingly enough, the following problem, apparently very closely related 
to the word problem, turns out to belong to the class of NP-complete prob- 
lems (recall that the Non-deterministic Polynomial class contains decision 
problems for which a guess solution can be checked in polynomial time; a 
particular NP problem is 'complete' if every other problem in the class can 
be polynomially reduced to it). 

Problem 2A \ Given a word (3 in the standard generators cri,cj2, 
. . . , cr„_i and their inverses, determine whether there is a shorter 
word (3' which represents the same element in B„. 




(8) 
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Finally 



2B. The best known algorithm for the conjugacy problem is ex- 
ponential in both \/3\ and n. 

It is worth noticing that the difficulty of solving the conjugacy problem 
in braid groups has been exploited for the construction of a public-key 
cryptosystem |TB] . 

3 Polynomial invariants of knots via braids 

Invariants of knots (links) of polynomial type arise (or can be reformulated) 
by resorting to representations of the braid group. Generally speaking, in 
order to represent the finitely presented group B„ defined in (jlj), we need 
an 'algebra' structure A, namely a vector space over some field (or ring) A, 
endowed with a multiplication satisfying associative and distributive laws. 
The algebra must have a unit with respect to multiplication and for our 
purposes must be also finitely generated, namely its elements can be decom- 
posed in terms of some finite 'basis' set, the number of elements of which 
equals the braid index n. 

The reason for considering an algebra should become clear if we recognize, 
on the one hand, that we can multiply braids € B„ by simply composing 
their diagrams: given Pi and P2 & B„ we get the product Pi P2 by placing 
the braid Pi above P2 and gluing the bottom free ends of Pi with the top 
ends of P2 (this opearation was implicitly assumed in (jlJ) and ©, see also 
Fig. 121). On the other hand, the operation associated with 'addition' of 
braids can be defined in terms of formal combinations of the type aPi + 6/32, 
for any Pi,P2 € B„ and a,b £ A (the field of scalars associated with the 
algebra A). 

With these premises, a representation ofBn inside the algebra A is a map 

PA :Bn ^ A (9) 

which satisfies 

Pa{PiP2) = Pa{Pi)pa{P2) V/3i,/32GB„, (10) 

namely is a group homomorphism from B„ to the multiplicative group 
G C A of the invertible elements of A (in particular: p/\{e) = 1, where e is 
the identity element of B„ and 1 denotes the unit of A; pa{P~^) = [pAiP)]~^ , 
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V/3). By using the standard generators of B„ defined in Q, it suffices to 
define the map Q on the elements of {cij} 

Pk{(Ti) = gi C k, (i = l,2,...n- 1), (11) 

and extend Unearly its action on products and sums of braids. Any pair 
of contiguous elements gi and ^j+i must satisfy the Yang-Baxter equation 
associated with the representation p/\, namely 

gi 9i+i 9i = gi+i 9i gi+i (12) 

while gigj = gjgi for \i - j\ > 1. 

Once defined the representation pA we may also introduce associated 
matrix representations of some fixed dimension by representing A over 
the algebra of (A^ x A^) matrices with entries in the field A 

A — > M(A,A^). (13) 

If we restrict the domain of the above map to the group G C A of in- 
vertible elements, the assignment ()13() can be rephrased as the choice an 
A^-dimensional vector space V over A, and thus we have the natural iso- 
morphism 

M{A,N) ^ GIa{V, N), (14) 

where GLa (V, N) is the general linear group of non-singular, A-linear maps 
V ^V. 

Loosely speaking, if we associate with a braid /? G B„ a matrix M{0) 
obtained by means of a representation (|13|) of dimension N = n, then /? can 
be characterized by a scalar, namely the trace of M{(3) (the 'character' of the 
representation in the group-theoretic language). Such traces are candidates 
to be interpreted as invariants of links presented as closed braids, cfr. © 
in the previous section. 

A trace function over the algebra A is formally defined as a linear function 
over A and, by extension, over a matrix representation algebra ()13() 

A — > M(A,Af) ^ A (15) 

satisfying the property 

Tr ( M(/3) M'{[3') ) = TV ( M'(/3') M(/3) ). (16) 

for any M{(3) ,M'{I3') which are the images under p/\ of two braids P,(3' 
G B„. It can be shown that Tr(M(/3)) is a link invariant since it does not 
change under Markov move of type i) (defined in section 1), namely 

Tr (M(/3)) = Tr (M'(/3')) if /? and /3'are conjugate. (17) 
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In other words, link invariants arising from these 'Markov traces' are regular 
isotopy invariants (we leave aside the issue of the invariance under Markov 
move of type ii)). 

The general algebraic setting outlined above is the framework under- 
lying the constructions of both the Jones and the HOMFLY ^7] link 
polynomials. In particular: 

• the Jones polynomial of a link L, J{L;t), is the Markov trace of the 
representation of B„ inside the Temperley-Lieb algebra TLn{t) (|13j. 
section 2, |18|^. It is a Laurent polynomial in one formal variable t 
with coefficients in Z, namely it takes values in the ring A = t~^]] 

• the HOMFLY polynomial P{L; t, z) is obtained as a one-parameter 
family of Markov traces (parametrized by z) of the representation of 
B„ inside the Hecke algebra Hn(t) (^SJ; section 4, Jl]). It is a Laurent 
polynomial in two formal variables with coefficients in Z, namely A = 
Z[t^\z^^. 

The relevance of such invariants is statistical mechanics and in quantum 
field theory is widely discussed in |^ [S] . 

As will be discussed in the following section, the problem of evaluating, 
or approximating, the Jones polynomial on a quantum computing machine 
has called many people's attention in the last few months. In this respect, 
the purely algebraic approach seems somehow lacking in selecting unitary 
representation of the braid group B„ in a natural way. For instance, the 
approaches proposed in |S] and jJUl provide two different types of Hilbert 
space structures and associated unitary representations of the braid group 
by resorting to clever, but 'ad hoc' constructions. 

The approach proposed by the authors of jHl El relies on the physical back- 
ground provided by 'unitary' topological quantum field theories, so that the 
problem of selecting a suitable unitary representation simply does not exist. 
On the other hand, a number of different mathematical tools can be called 
into play (all of which giving rise to the same link invariants) but unfor- 
tunately no one of them looks simpler than the purely algebraic method 
addressed above. In the rest of this section we shall provide a plain presen- 
tation of the quantum group approach leaving aside many technical details 
which can be found in ^Hj, pUj . 

Note preliminarly that, since we are looking for unitary representations 
(matrices) to be associated with link invariants, the running variable of 
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the polynomials (including in particular the Jones polynomial) has to be 
a unitary complex number c G C, with |c| = 1. The commonly adopted 
variable is a complex, r-th root of unity, namely 

q = exp(27ri/r) ,r G N, r > 1 (18) 

and the idea is that, by letting r grow, the polynomial can be evaluated in 
more and more points lying on the unit circle in C. The upgraded notation 
for the Jones polynomial is 

J(L;g) G nq,q-\ (19) 

The invariant of an oriented link L we are going to address is an extension 
of the Jones polynomial (fT^ . denoted by 

^{L;q;ji,j2,---,jM), (20) 

depending on the root of unity q introduced in (|18() and parametrized by 
labels {ji, j2, ■ ■ ■ , Jm} (the 'colors') to be assigned to each of the M link 
components {-Z>i}i=i,2,...,M- From the point of view of equivalence of links, 
J(L; q;ji,j2, ■ ■ ■ iJm) turns out to be a 'regular isotopy' invariant, but it can 
be shown that the quantity 

gl/2 _ q-1/2 J(^;9;Jl'i2, • • • , Jm), (21) 

where w{L) is the writhe of the link L defined in ©, is invariant under any 
ambient isotopy transformation. 

The colored polynomials H21|) reduce to Jones' ()19|) when all the colors 
ji, j2, . . . , jM are equal to a same j, with j = 1/2, but are genuine general- 
izations as far as they can distinguish knots with the same Jones polynomial 
j21j . Even more crucially, these invariants are 'universal', in the sense that 
they arise from a number of historically distinct approaches, ranging from 
R-matrix representations obtained with the quantum group method, mon- 
odromy representations of the braid group in 2D conformal field theories 
and the quasi tensor category approach by Drinfeld up to 2>D quantum 
Chern-Simons theory. All these models share the common feature of being 
'integrable', and integrability is reflected by the presence of Yang-Baxter- 
like equations, encoding the algebraic structure of braid groups in disguise. 

The basic objects which enters into the definition of the invariants (|2U|) 
are oriented and colored links and braids. Recall that a link is oriented if 
all its components {-t'j}i=i,2,...,M are endowed with an orientation. Since L 
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can be thought of as the closure /? of an open braid /? G B„ for some n 
(Alexander's theorem of section 1) each strand of (3 (/?) inherits naturally 
an orientation, depicted in figures by an arrow. 

The assignment of 'colorings' can be carried out in two different ways, by as- 
signing a color either to each oriented link component {m = 1,2,..., M) 
or to each strand [j (i = 1, 2, . . . , n) of the associated oriented braid (3 G B„. 
Of course the braid index n is (much) greater than the number of link com- 
ponents M, as can be easily recognized from the samples given in Fig. |l]and 
Fig. El In our paper [0] we adopted the latter option (see also [221 ) ^^sX gives 
rise to unitary, n-dimensional representations of B„ inside the representa- 
tion ring of the quantum group SU{2)q defined below. Here we are going to 
illustrate the first choice which is technically simpler and stresses the role of 
the so-called i?-matrix, namely the set of representations of the crossings 
in link diagrams as 'braiding operators' acting over the representation ring 
of SU{2)q. Note however that both the choices of the colorings provide the 
same colored link invariants, possibly up to an overall normalization factor. 

Coming at last to the point, the representation ring of SU{2)g, denoted 
by 9\{SU{2) q), can be introduced following in the footsteps of the con- 
struction of S'C/(2)-representation theory. According to our previous nota- 
tion, the ground ring (in which the link invariants will take their values) 
is A = Z[q^^] C C, with q = exp(27rz/r) as in H18() . The elements of 
9^(5C/(2)g) are complex Hilbert spaces, invariant under the action of the 
group (recall that a vector space V is invariant under the action of a group 
G G xV ^ V, namely transformed vectors keep on belonging to V; such 
spaces are referred to as invariant G-modules). As happens for SU{2), it 
can be shown that ^{SU{2)q), is spanned by finite-dimensional SU{2)q— 
modules {T^-^}- In the case of SU{2) the labels {j} (the spin quantum 
numbers from the quantum mechanical point of view) run over all integers 
and half-integers {0, ^, 1, |, . . .}, each is characterized by its dimension 
{2j + \) and is irreducible (namely cannot be decomposed into a direct sum 
of invariant subspaces of lower dimensions). 

In the g-deformed case it can be shown that the 5C/(2)g-modules {V^} are 
irreducible if and only if the labels {j} run over the finite set {0, ^, 1, |, . . . , r}. 
Each , spanned by (2j + 1) vectors, can be characterized by a scalar G A, 
the g-integer [2j + l]g, where [n]g = (g"/^ _ q~n/2yi^^i/2 _ ^1/2^ f^^. ^ ^ 

a positive integer. Thus, for each choice of the integer r, we have a distin- 
gushed family of irreducible representations (irreps) of SU {2)q 

dr = { },=o,...,r ; ^ [2j + 1], (22) 

which makes y{{SU{2)q) a finitely generated ring. As in the case of SU{2), 
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the ring structure is explicitated in terms of the direct sum © and tensor 
product (8) of irreps, namely 

W © &^{SU{2)q) ifj,k< r 

W © V'' e $K(SC/(2)g) iij + k<r, (23) 

where the ranges of the labels have to be suitably restricted with respect 
to the standard case. The analog of the Clebsch-Gordan series, giving the 
decomposition of the tensor product of two irreps into a (truncated) direct 
sum of irreps, reads 

min{ji+j2,r-ji~j2} 

yji ©yj2 = v^. (24) 

J = |jl-J2| 

Note however that the ring 9\{SU{2)q) is much richer than its 'classical' 
/S'C/(2)-counterpart because SU{2)q can be endowed with a quasitriangular 
Hopf algebra structure. This means that, besides the standard operators © 
and (8) we can also introduce a comultiplication A : SU (2)g — > SU{2)q ® 
SU{2)q, an antipode map A : SU{2)q — > SU{2)q, a counit e : SU{2)q C 
and a distinguished invertible element 

R G SU{2)q^SU{2)q, (25) 

called the R-matrix. We do not insist any further on the explicit definitions 
of A, ^ and e, and refer to |19j (section 1) for a quite readable account. 
The far-reaching role played by the i2-matrix becomes manifest when we 
define its action on the tensor product of a pair of irreducible SU{2)q- 
modules in 9\{SU{2)q). Denoting by R the operator associated to R, we 
have 

R : W © — > V'' W, (26) 

where, according to the values of the labels j, k have to be suitably 
restricted. These R-operators will be referred to as braiding operators asso- 
ciated with the i2--matrix ()25() . If we further extend the action of R to the 
ordered product of three irreps (E> V'' (g)V^ by defining 

R^.^ = R © Id : {V^ © V'') © T/' — > {V'' © V^) © 

Rfc, = Id © R : © (y^' © y') — > © (y' © V^) , (27) 
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where Id is the identity operator on the corresponding factor, then it can be 
shown that these operators satisfy the quantum Yang-Baxter equation 

Rjfc Rfci Rjfc = Rfc/ Rjfc Rfc/ • (28) 

The adjective 'quantum' refers of course to the underlying quantum group 
setting, while it easily recognized that coincides with ©, the algebraic 
Yang-Baxter relation characterizing the braid group structure, if we perform 
the substitutions 

ordered triple (j k I) i — > ordered triple (1 2 3) 
R (braiding operator) i — > R (crossing) . 

The explicit expression of the braiding operator R (and of its inverse R^^) 
can be worked out explicitly by selecting orthonormal basis sets in the 
S'C/(2)g -modules , V^, for each admissible choice of the pair j, k. In such 
bases, all the braiding operators and are unitary. 

Having collected all the necessary algebraic ingredients, the colored in- 
variant (|20|) for an oriented link L with M components can be now consis- 
tently interpreted as a single, A-linear map 

J(L;g;ji,j2,...,jM) : 9^(5C/(2),) A, (30) 

where the choice of the integer r in the root of unity H18() is constrained 
by the requirement r > M, at least in the most general case (M distinct 
colors) . 

The prescription for working out J (L; q;ji,j2, ■ ■ ■ iJm) can be summa- 
rized as follows (compare also Fig. [SJ where the plat presentation of the 
Borromean link is depicted). 

• Present the link L = yj Li [i = 1,2,..., M) as the plat closure of a 
braid and choose an orientation for each component (depicted by an 
arrow). Assign to each component a (distinct) color 

Li ^ ji (i = 1,2,...,M). (31) 

• Insert two parallel horizontal lines Ai, A2 cutting the 'cap' and 'cup' 
portions of the diagram, respectively. This choice provide the diagram 
with an overall, downward orientation. 

The region of the diagram lying between Ai and A2 is an open braid 
whose strands inherit suitable labels from the colorings (|31j) . 
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Figure 5: A plat presentation of the oriented and colored Borromean link 
as a closed braid on six strands. The parallel straight lines 1 (= Ai) and 2 
(= A2) intersect the diagram in points to be associated with Hilbert spaces 
which inherit the 'coloring' from the corresponding strands. 

• Assign to the intersection point between a line (Ai or A2) and the 
string labelled by j the SU{2)q irreducible module belonging to 
the distinguished family defined in p2jl . 

The whole configurations of intersection points on Ai and A2, each or- 
dered from left to right, are to be associated with the S'C/(2)g -modules 
y^i and , respectively, where each of them is the ordered tensor 
product of the individual irreps. To give an explicit expression of these 
correspondences, consider the case of the 3-components Borromean 
link depicted in Fig. |21 where in particular 

Vx, = ® (^V'' ®V'' ^V^ ^V^ 

Vx^ = V'^^V''® ®V^® V^. (32) 

Note that Vx^ and Vx2 have the same dimension as Hilbert spaces over 
C, given by the product of the dimensions of the individual factors. 
The number of such factors, say 2iV, is the same for the two spaces 
and equals the number of strands of the braid, or even the number of 
'caps' ('cups') lying above the line Ai (below A2) divided by two. This 
feature derives of course from the topological presentation we adopted 
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for the link L, since the braid obtained from the plat closure of any 
link has an even number of strands. 

• Going on with the example, in the representation ring y{{SU{2)q ) 
there exists a well defined, unitary operator B{L;q;j,k,l) to be asso- 
ciated with the trasformation relating V\-^ and in the diagram of 
the Borromean link L in Fig. [S] 

B{L;q;j,k,l):Vx, ^ Vx„ (33) 

where Vx-^ and Vx2 are explicitly defined in H32() . The composite braid- 
ing operator B(L; q; j, k, I) can be decomposed into an ordered sequence 
of the 'elementary' unitary braiding operators R (and their inverses) 
introduced in (|26|) . suitably tensorized with identities. The sequence 
is uniquely determined by going trough the diagram from Ai to A2. 

• In the case of the Borromean link, the matrix elements of the braiding 
operator (|33|) evaluated on (the tensor product of) orthonormal basis 
vectors of the spaces , can be collected into a unitary (2 J + 
1) X (2 J + 1) matrix parametrized by the colors j, k, I, namely 

U, k, I) G UiA, 2J + 1) (a, /? = 1, 2, . . . , 2 J + 1), (34) 

where U{A, 2 J + 1) is the algebra of unitary matrices on the ground 
ring A = and (2J + 1) = (2j + l){2k + 1)(2/ + 1). 

Finally, the colored link invariant J{L;q;j,k,l) is obtained by taking 
the trace of the matrix ()34|) , formally 

S{L-q-j,k,l) = (TrB,^)(j,fc,0, (35) 

where the resulting quantity turns out to contain the colorings through 
the 'quantum weights' [2j + l]q, [2/c+l]^, [2/+l]g (refer to the appendix 
of jl2] for the explicit expression of 

The above list of prescriptions can be applied to any link diagram as 
well, the 'output' being the colored link polynomial of the associated link, 
namely the (colored) trace of a suitable unitary matrix. Such 'Markov traces' 
are exactly the objects that can be handled in a quantum computational 
framework, as will be illustrated in the next section. In the present context 
the original Jones polynomial ()19() can be easily recovered by choosing the 
fundamental (^-spin) 2-dimensional irreducible representation on each of 
the link components {ii}i=i,...,M- 
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Note in conclusion that the whole construction does not depend on the 
particular plat presentation we choose for the link under examination, since 
the quantum Yang-Baxter equation (^5]) ensures that the braiding operators 
associated with different plat presentations of a same link can be converted 
one into the other. In other words, the role of this operatorial identity is 
specular to the role played by Markov move of type i) at the purely topolog- 
ical level. As indicated in (|2H) . the complete (ambient isotopy) invariance of 
the colored polynomials, implemented topologically by both Markov moves 
i) and ii), can be restored by taking into account the writhe of the link. 

4 Quantum computation of link polynomials 

Having defined in the second part of the previous section the SU (2)q-colored 
link polynomials, let us focus now on the Jones invariant (|19|) . which is the 
simplest, ^-spin colored polynomial, on the one hand, and the prototype 
of invariants arising in a purely algebraic context, on the other (c/r. the 
discussion at the beginning of section 3). The reason why Jones' case is so 
crucial in the computational context is actually due to the fact that a 'sim- 
pler' link invariant, the Alexander-Conway polynomial, can be computed 
efficiently, while the problem of computing 2-variable polynomials -such as 
the HOMFLY invariant briefly addressed in section 3- is NP-hard (see ^3] 
for the definitions of the mentioned invariants and j23j for an account of com- 
putational questions). The issue of computational complexity of the Jones 
polynomial in classical information theory can be summarized as follows. 

Problem 3 How hard is it to determine the Jones polynomial 
of a link L ? 

A quite exhaustive answer has been provided in [JH], where the evaluation 
of the Jones polynomial of an alternating link L at a root of unity q is shown 
to be 77^:P-hard, namely computationally intractable in a very strong sense 
(see the definition of this class below). 

A number of remarks are in order. Recall first that 'alternating' links are 
special instances of links, the planar diagrams of which exhibit over and un- 
der crossings, alternatively. Thus, the evaluation of the invariant of generic, 
non-alternating links is at least as hard. Secondly, the computation becomes 
feasible when the argument q of the polynomial is a 2nd, 3rd, 4th, 6th root of 
unity (refer to ^23^ for details on this technical issue). Recall finally that the 
T^P complexity class can be defined as the class of enumeration problems in 
which the structures that must be counted are recognizable in polynomial 
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time. A problem vr in ^^P is said #P-complete if, for any other problem vr' 
in #P, vr' is polynomial-time reducible to vr; if a polynomial time algorithm 
were found for any such problem, it would follow that ^^P CP. A prob- 
lem is T^P-hard if some #P-complete problem is polynomial-time reducible 
to it. Instances of #P-complete problems are the counting of Hamiltonian 
paths in a graph and the most intractable problems arising in statistical me- 
chanics, such as the enumeration of configurations contributing to ground 
state partition functions. 

The intractability of Problem 3 relies on the fact that it is not possible 
to recognize in polynomial time all the equivalent configurations of a same 
link L, namely link diagrams {D{L), D'{L), D"{L), . . .} related to each other 
by (regular) isotopy. Coming back to the problems addressed in section 2, 
since any link can be presented efficiently as a closed braid (Problem lA), 
we are justified in switching our attention to closed braids and implementing 
regular isotopy of diagrams by means of Markov move of type i). However, 
the intractability of Problem IB (selecting the diagram with the minimum 
braid index) prevents us from selecting an 'optimal' representation of the 
isotopy class of diagrams that would provide, in turn, a unique standard 
configuration to be handled for computational purposes. Moreover, owing 
to the fact that Markov move i) is closely related to the conjugacy problem 
in the braid group (Problem 2B), the whole matter could be reformulated 
within the framework of the theory of finitely presented groups as well. 
Then the issue of the optimal presentation turns out to be related also with 
the NP-complete problem stated in 2A' of section 2 (the 'shorter word' 
problem). 

In the discussion above, the relevant quantities encoding the 'size' of a 
typical instance of the computational problem -a link diagram L presented 
as a closed braid on n strands, L = /5 as in @- are of course the number 
of crossings k and the braid index n. We might consider, instead of At, the 
length 1/3 1 of the open braid f3 associated with L, which equals the number 
of generators (and inverse generators) in the explicit expression of /? as a 
word in B„. Finally, also the argument q of the Jones polynomial H19() is a 
relevant parameter since, when the integer r in (|18|) becomes ^ 1, we would 
reach more and more points on the unit circle in the complex plane, thus 
giving more and more accurate evaluations of the invariant. 

The computational intractability of Problem 3 does not rules out by any 
means the possibility of 'approximating' efficiently Jones invariant. 
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Problem 4- How hard is it to approximate the Jones polynomial 
J(L, q) of a link L at a fixed root of unity q (q ^ 2nd, 3rd, 4th, 
6th root)? 

Loosely speaking, the approximation we are speaking about is a number Z 
such that, for any choice of a small 5 > the numerical value of J{L,q), 
when we substitue in its expression the given value of q, differs from Z by 
an amount ranging between —6 and +6. In a probabilistic setting (either 
classical or quantum) we require that the value Z can be accepted as an 
approximation of the polynomial if 



where we agree to present always the link as a closed braid and refer the 
reader to [0] and JHI for more accurate statements of . 
In the framework of classical complexity theory there do not exists algo- 
rithms to handle Problem 4 and thus, at least at the time being, this problem 
is to be considered as intractable. From the quantum computational side, it 
was Michael Preedman who first addressed the general problem of evaluat- 
ing quantities of topological nature arising in the context of 3-D-topological 
quantum field theories and associated 2L'-conformal field theories [2^ (the 
physical quantum systems to be simulated are typically anyonic systems). 
Freedman's 'quantum field computer' was especially designed to this goal, 
although it was later recognized that this model of computation is actually 
(polynomial time) reducible to the standard quantum circuit model. The 
formal statement of the answer to Problem 4 was given in i2| (see also |2S1) 
and can be summarized as follows. 

4- The approximation of the Jones polynomial of a link pre- 
sented as the closure of a braid at any fixed root of unity is 
BPQ-complete. Moreover, this problem is universal for quan- 
tum computation, namely is the 'prototype' of all problems effi- 
ciently solvable on a quantum computer. 

Recall that BQP is the computational complexity class of problems which 
can be solved in polynomial time by a quantum computer with a probability 
of success at least ^ for some fixed (bounded) error. In [7j it was proved 
that P-* = BQP, where P"' is defined as the class of languages accepted in 
polynomial time by a quantum Turing machine with an oracle for the lan- 
guage defined by Problem 4. This equality between computational classes 
implies that, if we find out an efficient quantum algorithm for Problem 4, 



Prob 




(36) 
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then the problem itself is complete for the class BQP, namely each problem 
in this last class can be efficiently reduced to a proper approximate eval- 
uation of the Jones polynomial of a link (see fo^' ^ detailed discussion 
on this issue). According to the above remarks, the goal of working out 
efficient quantum algorithms for Problem 4 (not explicitly given in 0) does 
represent a breakthrough in quantum information theory. 

In the rest of this section we are going to illustrate, without entering in 
much technical details, the efficient quantum algorithms recently proposed 
by three groups, Aharonov, Jones and Landau 'W\, Garnerone, Marzuoli and 
Rasetti jHlE], Wocjan and Yard [TOl. 

Generally speaking, the approaches proposed in El and in |H1 El dif- 
fer both in the theoretical background underlying the construction of link 
invariants, and in the model of quantum computation used to deal with 
calculations. 

The Jones polynomial is defined in ^1 by resorting to the purely al- 
gebraic framework outlined in the first part of section 3, and these authors 
focus on the search for a unitary representation of the braid group B„ in the 
Temperley-Lieb algebra |H1 or in the Hecke algebra llOj . However, as already 
pointed out in section 3, such representations are selected on the basis of 
purely formal criteria, and the deep connection between Jones polynomial 
and topological quantum field theory P is left aside. 

In our approach the universality of Jones polynomial (or, better, of the col- 
ored link invariants (|2()|) ) in so many physical contexts makes it possible, 
not only to select a 'natural' unitary representation of the braid group, but 
justifies also the use of anyone of the equivalent mathematical backgrounds. 
In this paper we have described the colored polynomials in the algebraic 
framework provided by the representation ring yi{SU{2) g) (and associated 
R-matrix) since it does not require any previous knowledge of quantum field 
theoretic notions, not so familiar to mathematicians and computer scientists, 
either to experimental physicists. Fully equivalent methods based on 3D- 
Chern-Simons theory and associated Wess-Zumino-Witten 2L'-conformal 
field theory could have been used as well jSjini) making our approach closer 
to the spirit of Freedman's original paper |2ll. Moreover, the SU{2)q- 
colored link invariants are more general than the Jones polynomial with 
respect to the detection of knots ^|, on the one hand, and can be related 
in a quite simple way to more general topological invariants characteriz- 
ing closed three-dimensional spaces, on the other |141 1191 OU] . The latter 
remark opens the intriguing possibility of addressing, on an effective compu- 
tational ground, basic questions in field theories whose dynamical variables 
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are geometric objects, including quantum gravity models. 

For what concerns the computational frameworks, while the papers 
I1U| rely on techniques based on the standard quantum circuit model, in 
[SI IS] we exploit the spin-network quantum simulator proposed in to set 
up quantum-automaton-like implementations of braiding operators, each 
elementary step being worked in one unit of the intrinsic discrete time of 
the automaton itself. We are currently completing the proof that each of 
the latter elementary unitary operations can be 'efficiently' implemented 
also with respect to the 'standard' quantum circuit model 57_. 

The results of [S] can be summarized as follows (see also the recent 
review |28j). The unitary representation of the braid group B„ in the 
Temper ley-Lieb algebra (c/r. the first part of section 1 above for the al- 
gebraic background) is obtained by resorting to an adaptation of the 'path 
model' proposed earlier by Jones himself. The key idea is that each gen- 
erator of the braid group is mapped by the path model representation to 
a unitary matrix which can be simulated efficiently by a quantum circuit. 
In particular, each crossing can be implemented by a quantum circuit us- 
ing polynomially many elementary quantum gates, the resulting unitary 
matrix being no longer 'local' (namely, it operates non trivially on more 
than just two qubits corresponding to the two strands concurring at the 
crossing point), but in any case efficiently implemented. An entire braid 
(whose plat closure represents a link L) can be applied efficiently by em- 
ploying a number of elementary gates polynomial in the braid index n and 
in the number of crossings k. The approximation the Jones polynomial is 
reduced to the approximation of the Markov trace of the unitary matrix 
associated with the braid, which is carried out by a standard quantum al- 
gorithmic technique (the Hadamard test). In the work ^U] more general 
closures of braids (not just plat closures) are addressed, and the algorithm 
is based on a local qubit implementation of the unitary Jones- Wenzl repre- 
sentations of the braid group. 

In both the latter papers, the value of the integer r in the root of unity q is 
arbitrarily chosen, but constant, while a polynomial estimate of the growth 
of the time complexity function with respect to this parameter too has been 
established in [S] by resorting to field theoretic arguments. 

As already pointed out, our approach is quite different from the previous 
ones in working out the unitary matrix to be associated with the plat clo- 
sure of a braid, while the approximation of its trace, giving the colored link 
polynomials as in the example ()35() . can be carried out by similar standard 
techniques j27j . To enter in some more details, the 'physical' background 
provided by the 3D quantum SU{2) Chern-Simons field theory plays here 
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a prominent role, because the computational scheme of the spin-network 
simulator is actually designed as a discretized conterpart of the topo- 
logical quantum computation setting proposed in [241 125j . Our framework 
is designed to deal with the 5'C/(2)g-colored link polynomials (IHU]) viewed as 
vacuum expectation values of composite 'Wilson loop' operators in Chern- 
Simons theory, on the one hand, and with unitary representations of the 
braid group, on the other. These expectation values, in turn, provide a 
bridge between the theory of formal languages and quantum computation, 
once more having as natural arena for discussion the spin-network environ- 
ment. We actually implement families of finite states (and discrete time)- 
quantum automata capable of accepting the language generated by the braid 
group, and whose transition amplitudes are colored polynomials. In other 
words, our results can be interpreted in terms of 'processing of words' - 
written in the alphabet given by the generators of the braid group in the 
given representation- on a quantum automaton in such a way that the expec- 
tation value associated with the internal automaton 'evolution' is exactly the 
required link polynomial (after the application of a suitable 'trace operation' 
and the approximation of the trace within some fixed range as required in 
H36() ). Our quantum automaton calculation provides the unitary matrix as- 
sociated with (the plat closure of) a colored link L on n strands in a number 
of steps which is bounded linearly in the number of crossings k of the link, 
on the one hand, and polynomially bounded in terms of the braid index n 
(< nlnn), on the other. 
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